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Abstract 



Dynamical vacuum energy or quintessence, a slowly varying and spatially 
inhomogeneous component of the energy density with negative pressure, is 
currently consistent with the observational data. One potential difficulty with 
the idea of quintessence is that couplings to ordinary matter should be strongly 
suppressed so as not to lead to observable time variations of the constants of 
nature. We further explore the possibility of an explicit coupling between the 
quintessence field and the curvature. Since such a scalar field gives rise to 
another gravity force of long range (•> H^ 1 ), the solar system experiments 
put a constraint on the non-minimal coupling: |£| < 10 -2 . 

PACS numbers: 98.80.Cq; 04.80.Cc; 95.35.+d 
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I. INTRODUCTION 



Recently a number of observations suggest that the Universe is dominated by an energy 
component with an effective negative pressure 0. One possibility for such a component is 
the cosmo logical constant. Another possibility is dynamical vacuum energy or quintessence, 
a slowly varying and spatially inhomogeneous component with negative pressure @-|7[]. 

We face two problems when we consider such a nonzero vacuum energy. The first is the 
fine-tuning problem related to the energy scale of the vacuum energy density ~ 10~ 47 GeV. 
The second is the coincident problem: why the vacuum energy is beginning to dominate 
presently. While these two problems are separated in quintessence, they are degenerate for 
the cosmological constant, and one has to introduce the cosmo logical constant of extremely 
small energy scale at the very beginning of the universe. 

As a solution of the coincidence problem, the notion of a tracker field is introduced in 
||. It is shown that a very wide range of initial conditions approach a common evolutionary 
track, so that the cosmology is insensitive to the initial conditions similar to inflation. Once 
one parameter relating to the energy scale of the vacuum energy is fixed, the present-day 
equation-of-state wq = Pq/pq is automatically determined: there is a Qq — wq relation || . 

Direct methods to verify the idea of quintessence are important. Proposed possibilities 
are the following: the direct reconstruction of the effective potential from the luminosity 
distance - redshift relation observed for type la supernovae || ; the detection of quintessence 
from the measurements of a rotation in the plane of polarization of radiation from distant 
radio sources [|Hj . The direct interaction of the quintessence field to ordinary matter, how- 



ever, is found to be strongly suppressed so as not to violate the equivalence principle and 
the constancy of the constants of nature |[1(| . 

The possibility of an explicit coupling between the scalar field and the curvature is not 
excluded theoretically. It is then natural to consider further the coupling of the quintessence 
field to the gravity itself. In this paper, we examine the cosmological consequence of the 
non-minimal coupling of the quintessence field to the gravity. Since such a scalar field gives 
rise to both the time variation of the gravitational constant and a gravity force of long range, 
such a coupling should be constrained by experiments. 



II. NON-MINIM ALLY COUPLED QUINTESSENCE 

The action we consider is 

+ S m , (2.1) 

where k 2 = SnGbare is the bare gravitational constant and S m denotes the action of matter. 
The effective gravitational "constant" is defined by k 2 ^ = k 2 (1 — £k 2 2 ) -1 . £ is the non- 
minimal coupling between the scalar field and the curvature. In our conventions, £ = 1/6 
corresponds to the conformal coupling. 

We assume that the universe is described by the flat homogeneous and isotropic universe 
model with the scale factor a. The time coordinate is so normalized that a = 1 at the 
present. The field equations are then 
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where Pb,Pb denotes the background energy density, pressure, respectively, and V = dV/dcj). 
We consider a potential of inverse-power as an example of the tracker field for £ = 



V{4>) = M 4 (cj)/M)- a . (2.6) 
For £ = 0, there exists the following scaling solution during the background dominated epoch 
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The equation-of-state Wq is 



a + 2 



(2.7) 
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Since we consider a potential whose present mass scale is extremely small (^ Hq ~ 10 _33 eV), 
the force mediated by the scalar field is of long range, and hence the usual solar system limit 
on £, likewise the Brans-Dicke parameter, does apply. The correspondence to the Brans- 



Dicke field &bd and the coupling function u>($bd) of scalar-tensor theories of gravity JTT 
is given by 



8tt(1 - £k 2 2 )/k 2 
1 - ^ 2 2 



4£ 2 K 



2 2 /A2 



4£(8tt - k 2 $ sd ) 



(2.11) 
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A. perturbative analysis 

To consider the effect of £ qualitatively, we consider the case of |£|k 2 2 <C 1. Then during 
the background dominated epoch, Eq. ( [2.2|) and Eq. ( |2.4|) are approximated to 

H 2 = ^-p B (2.13) 

+ 3H<j) + tK 2 (l-3wB)pB<f) + V' = 0, (2.14) 

where we have used Eq.( |2.3| ) to derive Eq.( |2.14D . It is recently established that the scaling 
solutions Eq. (|2.71 - |2.8j ) with the same power- index persist even if £ ^ and that the stability 
of them does not depend on £ fl2 . 
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To the lowest order in £, the corresponding Brans-Dicke parameter is given by 

1 - £ K 2 2 3 1 



Aa{a + 2)i 2 '- 



where we have used the relation that holds for the potential of inverse-power 
the present-day value of the scalar field: 

<p z 2 ^ a 

Up to 0(£), the time variation of the gravitational constant is given by 
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G 
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(2.15) 
to estimate 

(2.16) 



(2.17) 



Hence, for £ > the gravitational "constant" increases with time, while it decreases for 
£ < 0. Eq. fl2.17| ) also shows that \G/G\ is larger for larger a since the potential then 
becomes steeper and the scalar field rolls down the potential more rapidly. 



B. constraining £ 

We perform the numerical calculation to examine in detail the time variation of G and the 
deviation from general relativity induced by the non-minimal coupling of the quintessence 
field to the curvature. Initial condition is set at a = 10~ 14 . We vary the fraction of the 
energy density of the quintessence field relative to radiation from 10 -9 to 10 -30 . We also 
choose various initial and <fi. We confirmed the tracking behavior: convergence to a 
common evolutionary track @,12]. Below we show typical results for the potential Eq. ( |2.6|) 



with a = 4. We choose the following parameters: = K lffPAi/3H 2 \ = 0.3 and H = 
100/ikm/sec/Mpc with h = 0.6. 

There exist a lot of experimental limits on the time variation of G [13|. Radar ranging 
data to the Viking landers on Mars gives \G/ G\ = (2±4) x 10 _12 yr _1 |14]]. Lunar laser ranging 
experiments yield \G/G\ = (0±11) x 10- 12 yr _1 and recently updated as \G/G\ = (l±8)x 
10~ 12 yr _1 ||16|| . More recently, a tighter bound is found by analysing the measurements of the 
masses of young and old neutron stars in binary pulsars: \G/G\ = (0.6±2.0) x 10~ 12 yr _1 |17| , 



although the uncertainties in the age estimation may weaken the constraint. Considering 
these experimental results, we will adopt the limit: \G/G\ = (0 ± 8) x 10 _12 yr _1 , and the 
limit by Thorsett is treated separately. 

In Fig. 1, we show the numerical results of G/G. The shaded region is already excluded 
by the current experimental limits. To examine the model dependencies of the results, we 
also show G/G for the potential of the form M 4 (exp(l//t</>) — 1) [Q by a dotted curve. We 
find that negative £ is severely constrained, while positive £ is loosely constrained and the 
limit is dependent on the potential. 

These results are intuitively understood via conformally transformed picture |18|]. If we 



perform the conformal transformation so that the scalar field is minimally coupled: 

^o6=^ii-«wr l - (2-i8) 



4 



Then the action Eq.( |2.1D becomes 



S = I d xy—g 

where 



^-^ 2 (0)(V0) 2 -^(0) 



+ S m , (2.19) 



1 - £ (1 - 6£)k 2 </> 2 
F ^ = » ( 2 - 2 0) 



V{<t>) 



Hence, after redefining the scalar field so that the kinetic term is canonical: 

$ = J d<f>F(cf>), (2.22) 

the action is reduced to that of the scalar field minimally coupled to the Einstein gravity. 
We can follow the dynamics qualitatively by simply looking at the effective potential V($>). 
Note that 1/(1 — £k 2 2 ) 2 is a decreasing function of <j) for £ < 0, while an increasing function 
for £ > 0. For £ < the effective potential V(<p) decreases more rapidly than V(0) (in 
particular, V($) decreases exponentially for large K<p), and consequently the scalar field 
rolls down the potential more rapidly. On the other hand, for £ > 0, V($) diverges at 
K 2 (j) 2 = so the slope of the effective potential becomes gentler and the scalar field rolls 
down the potential more slowly, and hence \G/G\ becomes smaller than that for £ < 0. 

We may summarize that the current experimental limits on the time variation of G 
constrain the non-minimal coupling as 

- 1(T 2 <£ < 1(T 2 ~ HT\ (2.23) 

while if the tighter limit by Thorsett is adopted, then we have 

- 1(T 2 < £ < 1(T 2 . (2.24) 

However, the limit is sensitive to the shape of the potential. 

Most important experimental limits we must consider are the solar system experiments, 



such as the Shapiro time delay and the deflection of light [19| because the non-minimally 
coupled scalar field can mediate the long range gravity force in addition to that medi- 
ated by a metric tensor. The recent experiments set constraint on the parameterized-post- 
Newtonian(PPN) parameter 7ppn as 



|7ppn - 1| < 2 x 10" 3 , (2.25) 
which constrains the Brans-Dicke parameter through the relation 7ppn = (ou + l)/(u + 2)|o 

El 

ujq > 500. (2.26) 

In Fig. 2, we show the present-day Brans-Dicke parameter defined by Eq. (|2.12| ) as a 
function of £. We also plot a curve derived under the assumption of |£|k 2 2 1, Eq. fl2.15| ). 



5 



We find a good agreement. Thus, using Eq. ( |2.15|) and Eq.( |2.26|) , the non-minimal coupling 
f is found to be constrained as 



|f | < 3.9 x 10" 2 = < 2.2 x 1(T 2 (2.27) 

fa(a + 2) 



as long as a > 1. The limit is less sensitive to the potential than that derived from \G/G\ 
because u does not explicitly depend on <fi unlike \G/G\. We note that the limit is found 
to be insensitive to Qm as long as VLm ~ 0.7. There is another PPN parameter /?ppn which 
is written in terms of to as (3ppN — 1 = (doj / d<& BD ){2oj + 4) _1 (2u; + 3)~ 2 | jl9| . The most 
recent results of the lunar- laser-ranging [^TJ, combined with Eq. (|2.25|) , yields 

|/?ppn - 1| < 6 x 10" 4 . (2.28) 

We find that |/?ppn — 1| ~ C(f 3 ) and consequently the experimental limit on Pppn is always 
satisfied if the condition Eq. (|2.25|) is satisfied. 



III. SUMMARY 



We have explored the possibility of an explicit coupling between the quintessence field 
and the curvature. Because the force mediated by the scalar field is of long range(> H^ 1 ), 
such a coupling is constrained by the solar system experiments. Through both analytical 
estimate and numerical integration of the equations, we have found that the limit is given 
by If I ~ 10~ 2 . The current limit on the non-minimal coupling, |f | < 10~ 2 , is not so strong 
when compared with other couplings to ordinary matter. For example, a coupling with 
the electromagnetic field is suppressed at the level of ^ 1CT 6 ; the coupling with QCD is 
at most ^ 10~ 4 [[Tod . We have also found that the induced time variation of G is sensitive 
to the shape of the potential. The future improvements in the limit of G/G may further 
constrain negative f or might lead to a detection of G/G < depending on the potential of 
the quintessence field. 
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FIG. 1. The present-day G/G as a function of £ for the potential of inverse power with 
a = 4(solid curve) and for the exponential potential(dashed curve). An approximated relation for 
|£|k 2 </> 2 -C 1 (Eq.( 2.17 )) is plotted as a dotted line. The shaded region is already excluded by the 
current experimental limits: |G/G| = (0 ± 8) x 10 -12 yr -1 . The limit by Thorsett is also shown: 



\G/G\ = (0±2) x l(T 12 yr -1 [|17| 
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Fig.2 

FIG. 2. The present-day Brans-Dicke parameter. The limit by the solar system experiments is 
ujq > 500 The solid curve is for the potential of inverse power with a = 4; the dashed curve 
is for the exponential potential. An approximated relation for |£|av 2 </> 2 <C 1 (Eq.( |2.15| )) is plotted 
as a dotted curve. 
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